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Abstract. Thanks to the new generation of gravitational wave detectors LIGO and VIRGO, the 
theory of general relativity will face new and important confrontations to observational data with 

■ unprecedented precision. Indeed the detection and analysis of the gravitational waves from compact 
binary star systems requires beforehand a very precise solution of the two-body problem within 

. general relativity. The approximation currently used to solve this problem is the post-Newtonian 

one, and must be pushed to high order in order to describe with sufficient accuracy (given the 
sensitivity of the detectors) the inspiral phase of compact bodies, which immediately precedes their 
final merger. The resulting post-Newtonian "templates" are currently known to 3.5PN order, and are 
used for searching and deciphering the gravitational wave signals in VIRGO and LIGO. 
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O 1. INTRODUCTION 

cn 

o 

\Q I A compelling motivation for accurate computations of the gravitational radiation field 

O ■ generated by compact binary systems (i.e., made of neutron stars and/or black holes) is 

O , ■ the need for accurate templates to be used in the data analysis of the current and future 

generations of laser interferometric gravitational wave detectors. It is indeed recognized 
that the inspiral phase of the coalescence of two compact objects represents an extremely 
important source for the ground-based detectors such as LIGO and VIRGO, provided 
that their total mass does not exceed say 10 or 20 Mq (this includes the very interesting 
^ I case of double neutron-star systems), and for space-based detectors like LISA, in the 

■ case of the coalescence of two galactic black holes, if the masses are within the range 
between say 10^ and 10^ M©. 

For these sources the post-Newtonian (PN) approximation scheme has proved to be 
the appropriate theoretical tool in order to construct the necessary templates. A program 
started long ago with the goal of obtaining these templates with 3PN and even 3.5PN 
accuracy.^ Several studies, e.g. fE 111], have shown that such a high PN precision is 
probably sufficient, not only for detecting the signals in LIGO/VIRGO, but also for 
analyzing them and accurately measuring the parameters of the binary (such high- 
accuracy templates will also be of great value for detecting massive black-hole mergers 
in LISA). The templates have been first completed through 2PN order [3]. The 3.5PN 



To appear in the Proceedings of the Spanish Relativity Meeting "A Century of Relativity Physics" 
(ERE05), Edited by Lysiane Mornas and Joaquin Diaz-Alonso. 

^ Following the standard custom we use the qualifier nPN for a term in the wave form or (for instance) the 
energy flux which is of the order of 1 /c^" relatively to the lowest-order Newtonian quadrupolar radiation. 



accuracy (in the case where the compact objects have negligible intrinsic spins) has been 
achieved more recently [4*^. 

The calculation of the 3PN order turned out to be very intricate and quite subtle. 
The first step has been to compute all the terms, in both the 3PN equations of mo- 
tion jmislSD^ and 3.5PN gravitational radiation field lulIiKil, by means of the 
Hadamard self- field regularization 1 14. il5i] . A regularization is needed in this problem 
in order to remove the infinite self-field of point masses. However, a few terms were left 
undetermined by Hadamard's regularization, which correspond to some incompleteness 
of this regularization occurring at the 3PN order. These terms could be parametrized 
by some unknown numerical coefficients called ambiguity parameters. The second step 
has been to use the more powerful dimensional regularization 10] > which is technically 
based on analytic continuation in the dimension of space, which finally enabled to fix 
the values of all the ambiguity parameters lUTL flB*. 5, 19]. 

In Section El of this article we review and comment on the striking appearance of 
Hadamard self-field regularization parameters at 3PN order, and on their computation 
using dimensional regularization. Section |3] is devoted to the notion of the multipole 
moments of an isolated post-Newtonian extended source, at the basis of the construction 
of gravitational-wave post-Newtonian templates. In Section |4l we present two checks 
of the values of the latter ambiguity parameters, coming from the comparison between 
the binary's dipole moment and its center-of-mass vector on the one hand, and based 
on an argument from classical field-theory diagrams on the other hand. Finally, in 
Section |5l we consider the limiting case where one of the masses is exactly zero, 
and the remaining one moves with uniform velocity, and show that such "boosted 
Schwarzschild solution" limit yields the determination of the third ambiguity parameter 
in the radiation field. These tests, altogether, provide a verification, independent of 
dimensional regularization, for all the ambiguity parameters in the 3PN gravitational 
radiation field. 

2. HADAMARD REGULARIZATION PARAMETERS 

The standard Hadamard regularization yields some ambiguous results for the computa- 
tion of certain integrals at the 3PN order, as Jaranowski and Schafer |0, 0] first noticed 
in their computation of the equations of motion of point particles within the ADM- 
Hamiltonian formulation of general relativity. Hadamard's regularization is based on 
the notion of partiefinie of a singular function, given by the angular integral of the finite 
part coefficient in the singular expansion of that function near a singular point, and the 
related notion of partie finie of a divergent integral. It was shown [6, 7] that there are 
two and only two types of ambiguous terms in the 3PN Hamiltonian, which were then 
parametrized by two unknown numerical coefficients ©static and COidnetic- 

Motivated by the previous result, Blanchet and Faye introduced an extended version 
of Hadamard's regularization 120. 21.1 . which is mathematically well-defined and free 
of ambiguities; in particular it yields unique results for the computation of any of 
the integrals occuring in the 3PN equations of motion. Unfortunately, the extended 
Hadamard regularization turned out to be in a sense incomplete, because it was found llsl 
1^ that the 3PN equations of motion involve one and only one unknown numerical 



constant, called A, which cannot be detemiined within the method. The comparison 
with the work [6, 7], on the basis of the computation of the invariant energy of compact 
binaries moving on circular orbits, revealed |8] that 



«stauc = -y^-^- (2) 

Therefore, the ambiguity COkinetic is fixed, while A is equivalent to the other ambiguity 
ecstatic- Notice that the value (d) for the kinetic ambiguity parameter (Ujiinetio which is in 
factor of some velocity dependent terms, is the only one for which the 3PN equations of 
motion are Poincare invariant. Fixing up this value was possible because the extended 
Hadamard regularization 1I20I I21I1 was defined in such a way that it keeps the Poincare 
invariance. 

The appearance of one and only one physical unknown coefficient A in the equations 
of motion constitutes a quite striking fact, that is related specifically with the use of some 
Hadamard-type regularization. Technically speaking, the presence of the parameter A is 
associated with the so-called "non-distributivity" of Hadamard's regularization. ^ Math- 
ematically speaking, A is probably related to the fact that it is impossible to construct 
a distributional derivative operator satisfying the Leibniz rule for the derivation of the 
product. The Einstein field equations can be written into many different forms, by shift- 
ing the derivatives and operating some terms by parts with the help of the Leibniz rule. 
All these forms are equivalent in the case of regular sources, but since the distributional 
derivative operator violates the Leibniz rule they become inequivalent for point particles. 
Finally, physically speaking, we can argue that A has its root in the fact that, in a com- 
plete computation of the equations of motion valid for two regular extended weakly self- 
gravitating bodies, many non-linear integrals, when taken individually, start depending, 
from the 3PN order, on the internal structure of the bodies, even in the "compact-body" 
limit where the radii tend to zero. However, when considering the full equations of mo- 
tion, we expect that all the terms depending on the internal structure can be removed, in 
the compact-body limit, by a coordinate transformation (or by some appropriate shifts 
of the central world lines of the bodies), and that finally A is given by a pure number, 
for instance a rational fraction, independent of the details of the internal structure of the 
compact bodies. From this argument (which could be justified by invoking the effacing 
principle in general relativity [1220) the value of A is necessarily the one we shall obtain 
below, Eq. ©, and will be valid for any compact objects, for instance black holes. 

The ambiguity parameter Wstatic^ which is in factor of some static, velocity- 
independent term, was computed by Damour, Jaranowski and Schafer lil7l] by means 
of dimensional regularization, instead of some Hadamard-type one, within the ADM- 
Hamiltonian formalism. Their result is 

«static = 0. (3) 



^ By non-distributivity we mean that the Hadamard regularization of a product of functions differs in 
general from the product of regularizations. 



As Damour et al. aim argue, clearing up the static ambiguity is made possible by the 
fact that dimensional regularization, contrary to Hadamard's regularization, respects all 
the basic properties of the algebraic and differential calculus of ordinary functions: 
associativity, commutativity and distributivity of point- wise addition and multiplication, 
Leibniz's rule, and the Schwarz lemma. In this respect, dimensional regularization is 
certainly better than Hadamard's one, which does not respect the distributivity of the 
product and unavoidably violates at some stage the Leibniz rule for the differentiation 
of a product. 

The ambiguity parameter A is fixed from the result Q and the necessary link Q 
provided by the equivalence between the harmonic-coordinates and ADM-Hamiltonian 
formalisms. However, A was also computed directly by Blanchet, Damour and Esposito- 
Farese |18] applying dimensional regularization to the 3PN equations of motion in 
harmonic coordinates (in the line of Refs. The end result, 

A = -i^, (4) 
3080' ^ ^ 

is in full agreement with Eq. Q. Besides the independent confirmation of the value of 
COstatic or A, the work [18] provides also a confirmation of the consistency of dimensional 
regularization, because the explicit calculations are entirely different from the ones of 
Ref. iITtII : harmonic coordinates are used instead of ADM-type ones, the work is at the 
level of the equations of motion instead of the Hamiltonian, a different form of Einstein's 
field equations is solved by a different iteration scheme. 

Let us comment here that the use of a self-field regularization, be it dimensional 
or based on Hadamard's partie finie, signals a somewhat unsatisfactory situation on 
the physical point of view, because, ideally, we would like to perform a complete 
calculation valid for extended bodies, taking into account the details of the internal 
structure of the bodies (energy density, pressure, internal velocity field). By considering 
the limit where the radii of the objects tend to zero, one should recover the same 
result as obtained by means of the point-mass regularization. This would demonstrate 
the suitability of the regularization. This program was undertaken at the 2PN order by 
Grishchuk and Kopeikin |23, 24] who derived the equations of motion of two extended 
fluid balls, and obtained equations of motion depending only on the two masses mi 
and m2 of the compact bodies. At the 3PN order we expect that the extended-body 
program should give the value of the regularization parameter A (maybe after some 
gauge transformation to remove the terms depending on the internal structure). Ideally, 
its value should be confirmed by independent and more physical methods. One such 
method is the one of Itoh and Futamase ll25L l2^ . who derived the 3PN equations of 
motion in harmonic coordinates by means of a particular variant of the famous "surface- 
integral" method introduced long ago by Einstein, Infeld and Hoffmann [27]. This 
approach is interesting because it is based on the physical notion of extended compact 
bodies in general relativity, and is free of the problems of ambiguities due to Hadamard's 
self-field regularization. The end result of Refs. [25, 26] is in agreement with the 
complete 3PN equations of motion in harmonic coordinates \^ ^ and, moreover, is 
unambiguous, as it does determine the ambiguity parameter A to exactly the value Q. 

We next consider the problem of the binary's radiation field, where the same phe- 
nomenon occurs, with the appearance of some Hadamard regularization ambiguity pa- 



rameters at 3PN order. More precisely, Blanchet, Iyer and Joguet lll2ll . in their compu- 
tation of the 3PN compact binary's mass quadrupole moment lij, found it necessary to 
introduce three Hadamard regularization constants ^, K and ^, which are additional to 
(and independent of) the equation-of-motion related constant A . The total gravitational- 
wave flux at 3PN order, in the case of circular orbits, was found to depend on a single 
combination of the latter constants, 6 = ^ + 2k+ i^, and the binary's orbital phase, for 
circular orbits, involves only the linear combination of and X given by = — 7A/3, 
as shown in 

Dimensional regularization (instead of Hadamard's) was applied in Refs. |H, [I^ to 
the computation of the 3PN radiation field of compact binaries, finally leading to the 
following unique values for the ambiguity parameters 

I = (5) 
5 9240 ' ^ ' 

K = 0, (6) 

C = -ry P) 

These values represent the end result of dimensional regularization. However, we shall 
review in the present Article some alternative calculations which provide some checks, 
independent of dimensional regularization, for all the parameters ©-O. 

The result ©-(jVl) completes the problem of the general relativistic prediction for the 
templates of inspiralling compact binaries up to 3PN order (and actually up to 3.5PN 
order as the corresponding tail terms have already been determined iTT]). The relevant 
combination of the parameters entering the 3PN energy flux in the case of circular orbits 
is now fixed to be 

The orbital phase of compact binaries, in the adiabatic inspiral regime (i.e., evolving 
by radiation reaction), involves at 3PN order a combination of parameters which is 
determined as 

7 , 1039 

e = e--?i = — — . (9) 

3 4620 

The fact that the numerical value of this parameter is quite small, 9 ~ 0.22489, indicates 
that the 3PN (or, even better, 3.5PN) order should provide an excellent approximation 
for both the on-line search and the subsequent off-line analysis of gravitational wave 
signals from inspiralling compact binaries in the LIGO and VIRGO detectors. 



3. THE MULTIPOLAR POST-NEWTONIAN FORMALISM 

3.1. Multipole moments of a post-Newtonian extended source 

The multipole moments of a post-Newtonian (PN) source, by which we mean a source 
which is at once slowly moving, weakly stressed and weakly self-gravitating, are crucial 
for the present gravitational wave generation formalism. They are obtained in Ref. [l28l] 



as functionals of the PN expansion of the pseudo-stress energy tensor t'^^ of the matter 
and gravitational fields in the harmonic coordinate system. The pseudo-tensor t^^ has a 
non-compact support because of the contribution of the gravitational field which extends 
up to infinity from the source. Let us denote the formal PN expansion of the pseudo 
tensor by means of an overbar, so that T^^ = PN[t^^]. The two types of multipole 
moments of the gravitating source, mass-type It moments and current-type ones Jl, 
are then given by ^ 



c^ b=qj y [(] 



c{i+l){2l + 3) '"[f+i 
2(2£+l) 



' c2(£+l)(£ + 2)(2£ + 5f ^'^'^ [ml'} ' ^^^^ 

^^^^ -XL-Dac^'i':]. (11) 



c(£ + 2)(2£ + 3) ""^^"'[£+1] 

Since Eqs. (frob - lfTTT) are valid only in the sense of PN expansions, the operational 
meaning of the underscript [£] in (fTUt- ifTTT) is actually that of an infinite PN series, which 
is given by 



r d 
dt 



- T^^(x,0, (12) 

(2£+l)!! 
" i2k)Ui2i + 2k+l)W 

A basic feature of the expressions of the moments is that the integral formally extends 
over the whole support of the PN expansion of the stress-energy pseudo-tensor, t'^^, i.e. 
from r = |x| = up to infinity. Recall that a formal PN series such as t'^^ is physically 
meaningful only within the near-zone. Therefore the integrals (fTOb - lfTTT) physically refer 
to a result obtained from near- zone quantities only (in the formal limit where c — > +°°). 
However, it was found extremely useful in Ref. [2^] to mathematically extend the 
integrals up to r ^ +oo. This was made possible by the use of the prefactor r^, together 
with a process of analytic continuation in the complex B plane. ^ This shows up in 
Eqs. (frob - lfTTT) as the crucial Finite Part (FP) operation, when 5 — > 0, which technically 



Our notation is: L = /i • • • i( for a multi-index composed of £ multipolar indices ii, - • • , if, xl = Xi^ ■ ■ -x,^ 
for the product of £ spatial vectors x' = jc,; and xi = STF(x,j • • -x,^ ) for the symmetric-trace-free (STF) 
part of that product, also denoted by carets surrounding the indices, jc^^j = i^. 

^ The prefactor should in principle be adimensionalized as [r/roY where ro is a constant arbitrary 
scale, but here we set ro = 1 . 



allows one to uniquely define integrals which would otherwise be divergent at their upper 
boundary, r = |x| +00. See Ref. Ii28i1 for the proof and details. 



3.2. Surface-integral expressions of the multipole moments 

Let us next review the recent derivation of an alternative form of the PN source 
moments (fTOb - dTTTl in terms of two-dimensional surface integrals. Such a possibility of 
expressing the moments, for general £ and at any PN order, as some surface integrals is 
quite useful for practical purposes, as we shall show in the application we consider in 
Section|5l In keeping with the fact that the "volume integrals" Eqs. (fTOb - ifTTT) physically 
involve only near-zone quantities, the "surface integrals" into which we shall transform 
the moments \i and J/^ physically refer to an operation which extracts some coefficients 
in the "far near-zone" expansion of the gravitational field, i.e. in the expansion in 
increasing powers of 1/r of the PN-expanded near- zone metric. Technically, as our 
starting point (fTOb - dTTT) is made of integrals extended up to r ^ +00, our mathematical 
manipulations below will involve "surface terms" on arbitrary large spheres r = Si. All 
these manipulations will be mathematically well-defined because of the properties of 
complex analytic continuation in B. 

The basic idea is to go from the "source term", T^^, to the corresponding "solution" 
, via integrating by parts the Laplace operator present in the Einstein field equation 

in harmonic coordinates, namely f^"^ = -j^^D/z^^, where h^"^ is the (PN expansion 
of the) basic gravitational field variable, satisfying the harmonic-coordinate condition 
dyh^'' = 0. From Eq. O we have 

„ 4 -|-oo / 7 \2k „ 

in the right-hand-side of which we insert □ = A — j , and operate the Laplacian 

by parts using A(r^+2fc^^) ^ {B + lk){B + U + Ik + \)r^+^''^^ xl- In the process we 
can ignore the all-integrated surface terms because they are identically zero by complex 
analytic continuation, from the case where the real part of B is chosen to be a large 
enough negative number. Using the expression of the coefficients (fT3l) . we are next led 
to the alternative expression 

/4 +00 / T \2k„ 

d'xr^XL T^" = y B{B + 2i + 4k+ 1) a^, ( 4" I / d\r'+^'-^ XLh^\ 

[i] IbTtG ' \cdt J J 

(15) 

A remarkable feature of this result, which is the basis of our new expressions, is the 
presence of an explicit factor B in front of the integral. The factor means that the result 
depends only on the occurrence of poles, oc l/B^, in the boundary of the integral at 
infinity: r +00 with t = const. 

Thanks to the factor B we can replace the integration domain of Eq. (fT5l by some 
outer domain of the type r>^%, where ^ denotes some large arbitrary constant radius. 



The integral over the inner domain r < 1% h always zero in the limit 5^0 because the 
integrand is constructed from 'x^^ , and we are considering extended regular PN sources, 
without singularities. Now, in the outer (but still near- zone) domain we can replace the 

— U V 

PN metric coefficients h by the expansion in increasing powers of 1 /r of the PN- 
expanded metric, which is identical to the multipolar expansion of the PN-expanded 
metric, that we shall denote by ^(/j^^) . Hence we have 

J [I] l6TZGj^Q ' \cdt J Jr>PJ ^ ' 

(16) 

We want now to make use of a more explicit form of the far near-zone expansion 
^ (^^^) , whose general structure is known. It consists of terms proportional to arbitrary 
powers of 1 / r, and multiplied by powers of the logarithm of r; more precisely. 



a,b 



where a can take any positive or negative integer values, and h can be any positive 
integer: a G Z, G N. The coefficients (p^^ depend on the unit direction n = x/r and 
on the coordinate time t (in the harmonic coordinate system). The structure ( flTl) for the 
multipolar expansion of the near- zone (PN-expanded) metric is a consequence of the 
so-called matching equation 



J{{}f )=J{{h\'^), (18) 

— U V 

which says that the multipolar re-expansion of the PN metric h agrees, in the sense 
of formal series, with the near-zone re-expansion (also denoted with an overbar) of the 
external multipolar metric ^ {h^^) (see 128] for details). Inserting Eq. Wh into (fT6b . 
we are therefore led to the computation of the integral 

/ d^^r^^^^-^XLJ^{h''^)=Yl^'^drr^^^^^'-'^{\nri' [ d^nLcp'^M , (19) 

where dQ. is the solid angle element associated with the unit direction n (and hi = xi/r^). 
The radial integral can be trivially integrated by analytic continuation in B, with result 



+00 / /7 \ ^ r (:S>B+2k+l-a+l 

^,^+2k+e-a^^^^^b^_fd\ ^ 



dB 



B + 2k + £-a + l 



(20) 



Remember that we are ultimately interested only in the analytic continuation of such 
integrals down to B = 0. And as an integral such as (l20b is multiplied by a coefficient 
which is proportional to B, we must control the poles of Eq. (l20b at 5 = 0. Those poles 
are in general multiple because of the presence of powers of Inr in the expansion, and 
the consecutive multiple differentiation with respect to B shown in Eq. (OUt . The poles 
at 5 = clearly come from a single value of a, namely a = 2k + i+l. For that value, the 



"multiplicity" of the pole takes the value b+l. Here a useful simplification comes from 
the fact that the factor in front of the integrals in (fT^ is of the form ~ 5(5 + ^) . In other 
words, this factor contains only the first and second powers of B. Therefore, only the 
simple and double poles l/B and 1/5^ in (l20b can contribute to the final result. Hence, 
we conclude that it is enough to consider the values b = 0,1 for the exponent of In r in 
the expansion (flTb . 

To express the result in the most convenient manner let us introduce a special no- 
tation for some relevant combination of coefficients 9^^(11, which as we just said 
correspond exclusively to the values a = i + 2k+ \ and Z? = or 1 . Namely, 



CCkl 



(2£ + 4^+l)(P27+,+i.o(n,0 + <' 



+f+l,lV"' 



(21) 



in which we have absorbed the numerical coefficient ttyt/ defined by dT^ . With this 
notation we then obtain 



FP B(B + 2i + 4k+l)ake 

B=0 



2k-2 



(22) 



where the brackets refer to the spherical or angular average (at coordinate time t), i.e. 



it) 



(23) 



The quantities (E^ are integrals over a unit sphere, and can rightly be referred to as 
surface integrals. These surface integrals are the basic blocks entering our alternative 
expressions for the multipole moments. If we wish to physically think of them as 
integrals over some two-surface surrounding the source, we can roughly consider that 
this two-surface is located at a radius with a <^ ^ <^ cT. Anyway, the important 
point is that, as we can see from Eq. (l23t . the surface integrals, and therefore the 
multipole moments, are strictly independent of the choice of the intermediate scale ^ 
which entered our reasoning. 

Finally, we are in a position to write down the following final results for an alternative 
form of the source multipole moments (fTOb - dmi . expressed solely in terms of the surface 
integrals of the type (l23t . 



Il 



■i-i' I / d \ 



— E — 



2k 



k=o k \cdt J 



(£+l)(2£ + 3) \cdtj 
2{2£+V 



4) (n^xnVi 



Jl 
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AG 



{i+l){i + 2){2i + 5) \cdtj 



kI+2 



(24) 



(25) 



4. MULTIPOLE MOMENTS OF TWO-BODY SYSTEMS 



4.1. Quadrupole and dipole moments, and the center-of-mass vector 

Let us show how a particular combination of ambiguity parameters can be determined 
within Hadamard's regularization and confirm the result of dimensional regularization. 
For this purpose we use the computations in Ref. [13] of the mass-type quadrupole 
lij and dipole I, moments of point particle binaries at the 3PN order. These were 
derived by applying the expression (flOb [with i = 1,2] to a binary systems of point 
masses, following the rules of the Hadamard regularization, in the so-called "pure 
Hadamard-Schwartz" (pHS) variant of it. Following the definition of Ref. [18], the pHS 
regularization is a specific, minimal Hadamard-type regularization of integrals, based on 
the usual Hadamard partie finie of a divergent integral, together with a minimal treatment 
(supposed to be "distributive") of compact-support terms. The pHS regularization also 
assumes the use of standard Schwartz distributional derivatives [15]. 

We shall denote by I-^^^ the result of such pHS calculation of the mass-type 
quadrupole moment. Now it was argued in Ref. lll2ll that the Hadamard regularization 
of the 3PN quadrupole moment is incomplete, in the sense that the pHS calculation ifj^^ 
must be augmented, in order to be correct, by some unknown, ambiguous, contributions. 
The first source of ambiguity is the "kinetic" one, linked to the inability of the Hadamard 
regularization to ensure the global Poincare invariance of the formalism. As discussed 
in Ref. [12] (see also Section |2l) we must account for the kinetic ambiguity by adding 
"by hands" a specific ambiguity term, depending on a single ambiguity parameter called 
^. The second source of ambiguity is "static". It comes from the a priori unknown 
relation between some Hadamard regularization length scales, s\ and $2 (one for each 
particles), and the ones, called r'^dsid rj, parametrizing the final 3PN equations of 
motion in harmonic coordinates |H,|a]. The static ambiguity is accounted for by two 
other ambiguity parameters and K (see Section |2l). 

The Hadamard-regularized 3PN quadrupole moment reads 



H^^K^Q = ir (26) 



+ 



rPHS 
ij 



c6 



1 mi-fm2\ {i j) , f ^ , 9 \ {ij) 



+ 1 



where one sees in the second term the effect of adding the ambiguities, parametrized 
by the same parameters t,, k and ^ as introduced in Ref. lfl2ll . Here, mi and m2 are the 
masses, y\, v\ and a\ denote the position, velocity and Newtonian acceleration of the 

first particle, and we pose y^^a^ = STF(yjaj) and vj'^^ = STF(vjv{). The symbol 1^2 
refers to the same terms but concerning the second particles. All the terms composing the 
pHS part have been explicitly computed up to 3PN order for general binary orbits 11 1311 . 

Let us now consider the case of the mass dipole moment I,-. Repeating the same 
arguments as for the quadrupole, we can write 1/ as the pHS part 1^^^^ and augmented by 
an ambiguous part. However, in the dipole case we find that no ambiguity of the kinetic 
type occurs, and that the only ambiguity is static. We find that the expression analogous 



to (E^ reads 



1,11 + ^] = ,PH» + 1 ^ + ^+ _L) + 1 „ 2. (27) 

As we see, there is only one ambiguity parameter, in the form of the sum of ^ and fc, 
where and K are exactly the same as in the quadrupole moment (l26b . Let us now fix 
that particular sum of ambiguity parameters. 

The case of the dipole moment I, is very interesting. Indeed let us argue that I,, which 
represents the distribution of positions of particles as weighted by their gravitational 
masses mg, must be identical to the position of the center of mass G,- of the system 
of particles (per unit of total mass), because the center of mass G, represents in fact the 
same quantity as the dipole I,- but corresponding to the inertial masses m\ of the particles. 
The equality between mass dipole I, and center-of-mass position G, can thus be seen as 
a consequence of the equivalence principle mi = mg, which is surely incorporated in 
our model of point particles. Now the center of mass G,- is already known at the 3PN 
order for point particle binaries, as one of the conserved integrals of the 3PN motion in 
harmonic coordinates. ^ The point is that G/, given in Ref. [10], is free of ambiguities; 
for instance the ambiguity parameter A in the 3PN equations of motion disappears from 
the expression of G/. Let us therefore impose the equivalence between and G,, which 
means that we make the complete identification 

I,[<^ + k]=G,. (28) 

Comparing I, with the expression of G, given by Eq. (4.5) in flo'], we find that Eq. (l28l) 
is verified for all the terms if and only if the particular combination of ambiguity 
parameters t,-\-K takes the unique value 

This result, obtained within Hadamard's regularization, is nicely consistent with the 
result of dimensional regularization, see Eqs. ©-Q. It shows that, although as we 
have seen Hadamard's regularization is physically incomplete (at 3PN order), it can 
nevertheless be partially completed by invoking some external physical arguments — 
in the present case the equivalence between mass dipole and center-of-mass position. 
On the other hand, dimensional regularization is complete; it does not need to invoke 
any external physical argument in order to determine the value of all the ambiguity 
parameters. Nevertheless, it remains that the result ( l29b . based simply on a consistency 
argument between the 3PN equations of motion and the 3PN radiation field, does provide 
a verification of the consistency and completeness of dimensional regularization itself. 



^ We neglect the radiation-reaction term at 2.5PN order 



4.2. Diagrammatic representation of the multipole moments 



Let us describe the multipole moments in terms of classical field-theory diagrams, 
representing the non-linear interactions of classical general relativity (we refer to llsoll 
for definition and use of these diagrams). We represent the basic delta-function sources 
entering the matter stress-energy tensor — i.e., the matter part of the pseudo-tensor 
t'^^ of Section|3] — as two world-lines, and each (post-Minkowskian) propagator D^^ 
as a dotted line. The various non-linear potentials entering the gravitational part of T^^ 
can then be represented by drawing some dotted lines which start at the matter sources, 
join at some intermediate vertices, corresponding to some non-linear couplings, and end 
at the field point x. Finally, we can represent the inclusion of the multipolar factors, such 
as XL, by adding a circled cross ®. It is then understood that one integrates over the 
crossed vertex, i.e., the field point. 

Using such a representation, the multipole moments are given by the sum of many 
diagrams. We are now looking at "dangerously" diverging diagrams, which generate 
poles 1 /e in a dimensionally continued approach, with d = 3 + £ being the dimension 
of space. Examining the types of singular integrals corresponding to the possible dia- 
grams, we find 1 19] that the only dangerously diverging diagrams are those containing 
(at least) three propagator lines that can simultaneously shrink to zero size, as a subset 
of vertices coalesce together on one of the particle world-lines. But as there are, in the 
present problem dealing with the 3PN order, at most three source points, this means that 
the dangerously divergent diagrams are only those represented in Fig.[T](or their mirror 
image obtained by exchanging 1 2). 



12 12 
a b 



FIGURE 1. Dangerously divergent diagrams contributing to the 3PN multipole moments. The world- 
lines of particles 1 and 2 are represented by vertical solid lines, the propagator D^' by dotted lines, the 
source points by bullets, and the (g) symbol means a multiplication by a multipolar factor, such as xl, 
together with a spatial integration / t/''x. 

Since the dangerous divergencies associated with the vicinity of the first world-line 
(say) are entirely contained in the diagrams shown in Fig. [H they are, therefore, pro- 
portional to mj (i.e., one factor mi per source point), without any explicit' dependence 
on the second mass m2. As a consequence, we can prove LI 9.] . because the presence 



^ There is also an implicit dependence on m2 via the fact that the acceleration a\ is proportional to t7i2- 
But, at the level of the diagrams, a\ must be considered as a pure characteristic of the first world-line. 



of ambiguity parameters is directly linked with the occurence of poles °^ 1/e, that the 
structure of the ambiguous terms in the mass quadrupole moment (E^ must be such that 
it is proportional to some factor mj. Now, the definition of the parameter K in Ref. [12] 
was to parametrize a conceivable a priori static ambiguity appearing in the renormaliza- 
tion of the logarithmic divergencies of the quadrupole moment, and these ambiguities 
were found in Eq. (l26b to be of the form + K") + Km\m2 (for what concerns the 
first particle). This shows that the parameter K corresponds to a mixing between dia- 
grams with three legs on the first world-line (as in Fig.[l]) and diagrams having two legs 
on the first world-line and one on the second. Our diagrammatic study has shown that 
the latter diagrams have no dangerous divergencies, i.e., that they do not introduce any 
conceivable ambiguity. Therefore we conclude, confirming Eq. that 

K = 0. (30) 



5. FIELD GENERATED BY A SINGLE BODY 

As another application, making use of the explicit surface-integral formula and 
yielding another check of ambiguity parameters, we wish to compute the source-type 
multipole moments of a spherically symmetric extended body moving with uniform ve- 
locity. Remember that our formalism assumes, in principle, that we are dealing with 
regular, weakly self-gravitating bodies. We expect, because of the effacing properties of 
Einstein's theory [22], that our final physical results, especially when they are expressed 
as surface integrals like in can be applied to more general sources, such as neutron 
stars or black holes. Indeed, we are going to confirm this expectation in the simplest pos- 
sible case, that of an isolated spherically symmetric body which is known, by Birkhoff 's 
theorem, to generate a universal exterior gravitational field, given by the Schwarzschild 
solution. 



5.1. The boosted Schwarzschild solution 



Following Ref. 11291] we shall apply our formulas to a boosted Schwarzschild solu- 
tion (BSS). Actually, in order to justify our use of the ESS (in standard harmonic 
coordinates), we must dispose of a small technicality. This technicality concerns the 
non-uniqueness of harmonic coordinates for the Schwarzschild solution, even under the 
assumption of stationarity (in the rest frame) and spherical symmetry. Indeed, under 
these assumptions, and starting from the usual Schwarzschild radial coordinate, say rg, 
the (rest frame) radial coordinate of the most general harmonic coordinate system, say 
r = k{rs), must satisfy the differential equation (see, e.g., Weinberg [31], page 181) 



_d_ 
drs 



2GM \ dk 
c^ ^ I drs 



2k. 



(31) 



The standard solution of Eq. (OTT) . which is considered in all textbooks such as 113 ill , 
reads simply 

r = e''-^^\rs)=rs-^. (32) 



In the black hole case, the solution (1^ is the only one which is regular on the hori- 
zon, i.e. when rs = 2GM/c^. However, in the case of the external metric of an extended 
spherically symmetric body, regularity on the horizon is not a relevant issue. What is 
relevant is that the solution of the external problem (ISTT) be smoothly matched to a regu- 
lar solution of the corresponding internal problem. As usual, this matching determines 
a unique solution everywhere. In general, this unique, everywhere regular, solution will 
correspond, in the exterior of the body, to a particular case of the general, two-parameter 
solution of the second-order differential equation (1311) . The latter is of the form 

r = k^---^-\rs) = Ci l^rs - ^) ^C^^ri) , (33) 

where k^ij^ denotes the (uniquely defined) "radially decaying solution" of Eq. (IbTT) . and 
where C\ and C2 are two integration constants. Indeed, when considering the flat- space 
limit of Eq. dBTT) . it is easily seen that there are two independent solutions which behave, 
when r5 +0°, as r^ and r^^ respectively. An explicit expression for the decaying 
solution is ^ 

(2,2.4, ?^?^). (34) 

We can always normalize C\ to the value C\ = \. Then, with the above definitions, C2 
has the dimension of a length cubed. By considering in more detail the matching of 
the general solution of the harmonically relaxed Einstein equations at the 2PN level 
(see, e.g., the book by Fock [32], page 322), one easily finds that the second integration 
constant is of the order of C2 ~ {GM / c^Y a, where a denotes the radius of the extended 
body under consideration. It is also easily checked that the constant C2 parametrizes, at 
the linearized order, a gauge vector of the form (p' oc C2 <?/( 1 /r) , and can thus be referred 
to as a "gauge parameter". 

Contrarily to the multipole moments of stationary sources, which are geometric 
invariants (and can be expressed as surface integrals on a sphere at spatial infinity), 
the source multipole moments defined in Ref. [28] (and re-expressed in Section E] as 
surface integrals over spheres in some intermediate region, a ^ r ^ c T) are probably 
not geometric invariants. They are useful intermediate constructs, which allow one to 
compute physically invariant information, but their definition is linked to the choice 
of harmonic coordinates covering the source. There are also various gauge multipoles 
(denoted W^, X^, Y^, in Ref. [28]) which will influence, at some non-linear order, 
the values of the two sequences of physical multipoles: II, Jl- Therefore, one should 
expect that, at some non-linear order, the physical multipoles \i, Ji of a boosted general, 
harmonic-coordinate spherically symmetric metric will start to depend on the value of 
the gauge parameter C2. 



Here F{2,2,4,z) denotes a particular case of Gauss' hypergeometric function 

^, „ , , aj3 z a(a + I)j3(j3 + 1) 



Here, we are only interested in computing the quadrupole moment I,j of a boosted 
general spherically symmetric metric. We shall see below that the index structure of 
lij will be provided by the STF tensor product of the boost velocity V with itself, de- 
noted V^'V-''^ (assuming that the origin of the coordinates is at the initial position of 
the center of symmetry of the BSS). Therefore, any contribution to 1,^ coming from 
the gauge parameter C2 must contain, at least, the factors C2 and V^'V-'\ and also the 
total mass M. Taking into account the dimensionality of C2 ~ {GM / c^)^ a, which is 
that of a length cubed, it is easily seen that there is no way to generate such a contri- 
bution to lij. Therefore, we conclude that the source quadrupole moment of a boosted 
general, harmonic-coordinate spherically symmetric metric is strictly equal to the source 
quadrupole moment of a boosted standard harmonic-coordinate Schwarzschild solution, 
obtained by setting C2 = (and Ci = 1) in (l33l) . i.e. by choosing the standard harmonic 
radial coordinate 

In the following, we shall therefore consider only such a boosted Schwarzschild solu- 
tion (BSS) in standard form. We shall sometimes refer to the source of this solution as 
a black hole (though, strictly speaking, one should always have in mind some extended 
spherical star). For simplicity, we shall translate the origin of the coordinate system so 
that it is located at the initial position of the black hole at coordinate time t = 0. With 
this choice of origin of the coordinates all the current-type moments of the BSS are 
zero. We shall concentrate our attention on the mass-type quadrupole moment that 
we shall compute at the 3PN order. 

The BSS metric, written in terms of the Gothic metric deviation h^"^ = ^f—gg^"^ — 
7]^"^ , satisfying standard harmonic coordinates so d^h^"^ = 0, is best formulated in a 
manifestly Lorentz covariant way as 



where is the time-like unit four-velocity of the center of symmetry of the BSS, where 
is the space-like unit vector pointing from the BSS to the field point along the di- 
rection orthogonal (in a Minkowskian sense) to the world line of the BSS, and where 
r^ denotes the orthogonal distance to the world line (square root of the interval). How- 
ever, in our explicit calculations (done with the software Mathematica) it is preferable to 
employ a more "coordinate-rooted" formulation of the BSS metric, which is of course 
completely equivalent to (though less elegant than) Eq. (l35b . 

Le us denote by = (ct.x) the global reference frame, in which the black hole is 
moving, and by = (cT, X) the rest frame of the black hole — both and are 
assumed to be harmonic coordinates. Let x'{t) be the rectilinear and uniform trajectory 
of the center of symmetry of the BSS in the global coordinates x^, and V = (V) be the 
constant coordinate velocity of the BSS, 




v 



(35) 



dx'{t) 
dt 



(36) 



The rest frame is transformed from the global one x^^ by the Lorentz boost, 



(37) 



For simplicity we consider a pure Lorentz boost A^y{Y) without rotation of the spatial 
coordinates. As explained above, we can assume that the metric of the BSS in the rest 
frame takes the standard harmonic-coordinate Schwarzschild expression, which we 
write in terms of the Gothic metric deviation H^^, satisfying d^H^^ = 0. Hence, 
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(38) 

(39) 
(40) 



where M is the total mass, R = |X| and A^' = X^/R. A well-known feature of the 
Schwarzschild metric in harmonic coordinates is that the spatial Gothic metric H'j is 
made of a single quadratic-order term oc as shown in Eq. (l40b . The Gothic metric 
deviation transforms like a Lorentz tensor so the metric of the BSS in the global frame 
reads as 



h^''{x)=A%A\HP''{A-'x) 



(41) 



in which the rest-frame coordinates have been expressed by means of the global ones, 
i.e. X^{x) = (A^^)^yjc^, using the inverse Lorentz transformation. The only problem is 
to derive the explicit relations giving the rest-frame radial coordinate R and the unit 
direction A^' as functions of their global-frame counterparts r and n', of the global 
coordinate time t, and of the boost velocity V. For these relations we find 
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(42) 
(43) 



where 7 = (l — V^/c^) ''^ and (Vn) = V ■ n = V-'n-' is the usual Euclidean scalar 
product. The latter formulation of the BSS metric is well adapted to our calculations 
because we shall have to perform, for computing the source multipole moments, an 
integration over the coordinate three-dimensional spatial slice x G M^, with coordinate 
time t = const, which is easily done using the explicit relations (l42b - (l43l) . 



5.2. Quadrupole moment of a boosted Schwarzschild black hole 

We compute the quadrupole moment 1;^ of the BSS, following the prescriptions 
defined by Eq. (l24ll . To this end we first expand h^^ when c +0°, taking into account 



all the c's present both in the expression of the rest frame metric H^^ as well as 
those coming from the Lorentz transformation (BTt -lB^l. In this process the boost 
velocity V is to be considered as a constant, "spectator", vector. Note in passing that, 
in the present problem, the characteristic size a of the source at time t is given by the 
displacement from the origin, a ~ Vf, where V = |V|, while the near- zone corresponds 
to r -C cL Therefore, the far near- zone, where we read off the multipole moments as 
some combination of expansion coefficients is the domain V t <ti r <ti ct. We 

have evidently to assume that V <^c for this region to exist. 

— uv 

We then first get the near- zone (or PN) expansion of the BSS metric, h , by expand- 
ing in inverse powers of c up to 3PN order. Next we compute the multipolar (or far) 
re-expansion of each of the PN coefficients when r +0° with t = const. In this way 
we obtain what we have denoted by ^(/?^^) in Eq. (flTt . In the BBS case it is evident 
that the far-zone expansion given by ST% involves simply some powers of 1 /r, without 
any logarithm of r. 

With^(^^'') in hand we have the coefficients of the various powers of 1 /r, and we 
obtain thereby the needed quantities defined by Eq. (ETT) . It is then a simple matter 
to compute all the required angular averages present in the formula (|24l) and to obtain 
the following 3PN mass quadrupole moment of the BSS, 



I 



BSS 



ij 



+ 



4 G^M^ 



1 ^Yl ^Yl 507 

^ 14 168 1232 ~^ 



1 



(44) 



The first term represents the standard Newtonian expression, augmented here by a bunch 
of relativistic corrections. (Recall that we have chosen the origin of the coordinate 
system at the initial location of the BSS at f = 0.) 

The last term in Eq. (l44l) . with coefficient ^ = 4/7, is the most interesting for our 
purpose. It is purely of 3PN order, and it contains the first occurence of the gravitational 
constant G, which therefore arises in the quadrupole of the BSS only at 3PN order. 
This term is interesting because it corresponds to one of the regularization ambiguities, 
due to an incompleteness of Hadamard's self-field regularization, which appears in 
the calculation of the mass-type quadrupole moment of point particle binaries at the 
3PN order 1. 12. ,13,1 . As we see from Eq. (E^ the associated ambiguity parameter is ^, 
which represents in fact the analogue of the kinetic ambiguity parameter (Okinetic in the 
equations of motion, see Eq. It is now clear that C, can be determined from what 
we shall now call the BSS limit of a binary system, which consists of setting one of the 
masses of the binary to be exactly zero, say m2 = 0. 

We have obtained the BSS limit of the 3PN mass -type quad rupole moment of compact 
binaries computed for general binary orbits in Refs. lll2Lll3h . We have also inserted for 
the position of the first body y\ = v\ t in order to conform with our choice for the origin 
of the coordinates. In this way we get 
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(45) 



The comparison of Eqs. (B31) and (R^ reveals a complete match between the two results 
if and only if we have the expected agreement between the masses, i.e. M = m\, and the 
velocities, v\ = V (since the velocity of the body remaining after taking the BSS limit 
should exactly be the boost velocity), and the ambiguity constant ^ takes the unique 
value 



Our conclusion, therefore, is that the ambiguity parameter ^ is uniquely determined 
by the BSS limit. Because of the close relation between the BSS limit with Lorentz 
boosts, it is clear that ^ is linked to the Lorentz-Poincare invariance of the multipole 
moment formalism of Ref. [28] as applied to compact binary systems in fT2',n"3]. This 
link strongly suggests that the specific value (l46t represents the only one for which the 
expression of the 3PN quadrupole moment is compatible with the Poincare symmetry. 
In other words the present calculation indicates that the Poincare invariance should 
correctly be incorporated into the laws of transformation of the source-type multipole 
moments for general extended PN sources as given by Eqs. (fTOb - ifTTTl or (l24ll - (l25t . 

Let us finally emphasize that Eq. (l46t has been obtained here without using any 
regularization scheme for curing the divergencies associated with the self field of point 
particles. However, we find, very nicely, that the value for C, is in agreement with the 
one derived in the problem of point particles binaries at 3PN order by means of the 
dimensional self-field regularization, Eq. Q. This shows in particular that dimensional 
regularization is able to correctly keep track of the global Poincare invariance of the 
general relativistic description of isolated systems. 
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